Probing CPT Violation with CMB Polarization Measurements 
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The electrodynamics modified by the Chern-Simons term Cca ~ p^A^F'^'^ with a non- vanishing 
Pfj, violates the Charge-Parity-Time Reversal symmetry (CPT) and rotates the linear polarizations 
of the propagating Cosmic Microwave Background (CMB) photons. In this paper we measure the 
rotation angle Aa by performing a global analysis on the current CMB polarization measurements 
from the seven-year Wilkinson Microwave Anisotropy Probe (WMAP7), BOOMERanG 2003 (BOS), 
BICEP and QUaD using a Markov Chain Monte Carlo method. Neglecting the systematic errors of 
these experiments, we find that the results from WMAP7, BOS and BICEP all are consistent and 
their combination gives Aq — —2.33 ± 0.72 deg (68% C.L.), indicating a Scr detection of the CPT 
violation. The QUaD data alone gives Aq = 0.59±0.42 deg (68% C.L.) which has an opposite sign for 
the central value and smaller error bar compared to that obtained from WMAP7, BOS and BICEP. 
When combining all the polarization data together, we find Aa = —0.04 ± 0.35 deg (68% C.L.) 
which significantly improves the previous constraint on Aa and test the validity of the fundamental 
CPT symmetry at a higher level. 
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Introduction - The accumulating high precision obser- 
vational data of the CMB temperature and polarization 
spectra are not only crucial to determine the cosmolog- 
ical parameters [ij, but also make it possible to search 
for new physics beyond the standard model of particle 
physics. One striking example along this line is the test 
of the CPT symmetry. As a fundamental requirement of 
particle physics, the CPT symmetry has been proved to 
be exact and well tested by various laboratory experi- 
ments. However, the validity of this symmetry needs to 
be reevaluated in the context of cosmology. And in fact, 
there have been some theoretical studies indicating the 
possible break-down of the CPT symmetry at some level, 
and interestingly, the cosmological measurements of the 
CMB polarization facilitate the direct detection of the 
CPT violating signal 

To begin with, consider an effective Lagrangian of elec- 
trodynamics including a Chern-Simons term Q Ccs ^ 
PfiAi^F^'^, where is an external vector, and F'^'^ = 
{l/2)e^^'^'"^Fpa denotes the dual of the electromagnetic 
tensor. Note that this model violates the Lorentz and 
CPT symmetries if is non-vanishing. Also, this effec- 
tive Lagrangian is not generally gauge invariant, but its 
action is invariant if = dfj,p^. This equality holds in 
some cases, for example, is a constant in spacetime or 
the gradient of a scalar field in the quintessential baryo- 
/leptogenesis Q ; or the gradient of a function of the Ricci 
scalar in the gravitational baryo- /leptogenesis @ . 

This CPT violating interaction yields a rotation, quan- 
tified by Aa, of the polarization vector of the electro- 
magnetic waves traveling over a distance on the cosmo- 
logical scale, and this mechanism is dubbed the Cosmo- 
logical Birefringence (CB) The rotation angle Aa is 
given in term of p^ by Aa ~ J Pfidx'^ 0], and it has 



imprints on the CMB polarization data, namely, all the 
CMB two-point functions, except for the temperature- 
temperature auto correlation (TT), will be altered, and 
most importantly, the cosmological birefringence can in- 
duce non-zero TB and EB spectra, which is vanishing in 
the standard cosmological model. Denoting the rotated 
quantity with a prime, one has the following relations 
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Given the CMB polarization data and Eq.(IT]), one can 
constrain the rotation angle to test the CPT symmetry. 

In this work, we report the latest result on the 
measurement of the rotation angle using the most 
up-to-date CMB polarization data including WMAP7, 
BOOMERanG 2003, BICEP and QUaD. 

CMB Polarization Measurements - In our previous 
analysis Q , we measured the rotation angle using the po- 
larization data from WMAP7 [lo| and the BOOMERanG 
dated January 2003 Antarctic flight [TJ]. The WMAP5 
polarization data are composed of TE/TB/EE/BB/EB 
power spectra on large scales (2 < £ < 23) and TE/TB 
power spectra on small scales (24 < £ < 800), while 
the BOS experiment measures the small-scale polariza- 
tion power spectra in the range of 150 < £ < 1000. 

Recently, the Background Imaging of Cosmic Extra- 
galactic Polarization (BICEP) [12l | and QU Extragalac- 
tic Survey Telescope at DASI {QVaD) [l3| collaborations 
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FIG. 1: The binned TB and EB spectra measured by 
the small-scale CMB experiments of BOOMERanG (black 
squares), BICEP (red circles) and QUaD (blue triangles). The 
black solid curves show the theoretical prediction of a model 
with Aa = -2.62 deg. 
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FIG. 2: One-dimensional posterior distributions of the rota- 
tion angle derived from various data combinations. The dot- 
ted vertical line illustrates the unrotated case (Aa = 0) to 
guide eyes. 



released their high precision data of the CMB tempera- 
ture and polarization including the TB and EB power 
spectra. These two experiments, locating at the South 
Pole, arc the bolometric polarimcters designed to cap- 
ture the CMB information at two different frequency 
bands of lOOGHz and 150GHz, and on small scales - 
the released first two-year BICEP data are in the range 
of 21 < ^ < 335 [lli; whereas the QUaD team measures 
the polarization spectra at 164 < I < 2026, based on an 
analysis of the observation in the second and third season 

In Fig[T]wc show the binned TB and EB power spec- 
tra released by the BOOMERanG, BICEP and QUaD 
collaborations. Compared to the B03 data, one can see 
that the BICEP and QUaD data have apparently smaller 
errors, implying that adding these data to the previous 
analysis, e.g. Ref.0], might in principal narrow down the 
constraint on the rotation angle, which is essentially the 
aim of this work. 

Method - Given the aforementioned CMB polarization 
data, we make a global analysis to constrain the rotation 
angle Aa using a modified version of CosmoMC , a p ublicly 
available Markov Chain Monte Carlo engine jl6[. With- 
out loss of generality, we assume the purely adiabatic ini- 
tial conditions and a flat universe, and explore the param- 
eter space of P = (wfc,a;c,0s,T,ns,log[lO^''^s],r, Aa). 
Here, a;;, = and Wc = are the physical baryon 

and cold dark matter densities relative to the critical den- 
sity, respectively, O^ denotes the ratio of the sound hori- 
zon to the angular diameter distance at decoupling, r 
measures the optical depth to re-ionization. As and Ug 
characterize the amplitude and the spectral index of the 
primordial scalar power spectrum, respectively, r is the 
tensor to scalar ratio of the primordial spectrum, and 



TABLE I; Constraints on the rotation angle from various 
CMB data sets. The Mean values and 68% C.L. error bars 
are shown. 



Data 




Aa (deg) 


Reference 


WMAP5H 


-B03-I-BICEP 


-2.62 ± 0.87 


This work 


WMAP7H 


-B03-I-BICEP 


-2.33 ± 0.72 


This work 


BICEP 




-2.60 ± 1.02 


This work 


WMAP7 




-1.1 ± 1.3 


Ref.[15j 


WMAP5H 


-BOS 


-2.6 ± 1.9 


Ref.[9] 


WMAP5 




-1.7±2.1 


Ref.[l] 


WMAP3H 


-BOS 


-6.2 ± 3.8 


Ref.[20] 


WMAP3 




-2.5 ± 3.0 


Ref.[21j 


WMAP3H 


-BOS 


-6.0 ±4.0 


Ref.[r| 



we choose fcso = 0.05 Mpc^^ as the pivot scale of the 
primordial spectrum. Furthermore, in our analysis we 
include the CMB lensing effect, which also produces B 
modes from E modes IVf, when we calculate the theo- 
retical CMB power spectra. 

The rotation angle Aa is accumulated along the jour- 
ney of CMB photons, and the constraints on the rotation 
angle depends on the multipoles I [l3|. In Ref.[H, the 
WMAP5 group found that the rotation angle is mainly 
constrained from the high-.^ polarization data, and the 
polarization data at low multipoles do not affect the re- 
sult significantly. Therefore, in our analysis, we assume 
a constant rotation angle Aa at all multipoles. Fur- 
ther, we also impose a conservative flat prior on Aa as, 
-7r/2 < Aa < 7r/2. 

Numerical Results ~ We present our result derived from 
the WMAP7, BOB and BICEP polarization data in Ta- 
bleUand FigHJ in comparison with the published results. 
Since it is still not very clear how to combine the system- 
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atic errors from different polarization measurements, in 
our calculations we do not include the systematic errors 
of the CMB measurements [l2l.[ll[l9j. 

As shown, the previously published constraints on the 
rotation angle, including the most stringent constraints 
Aa = -2.6 ± 1.9 deg (Icr) from WMAP5+B03 pre- 
sented in Ref.Q, are all consistent with Aa = at 
95% confidence level. However, in this work we find 
that the BICEP data alone give almost the same cen- 
tral value as that from WMAP5-I-B03, but tighten the 
constraints by roughly a factor of two, giving Aa 
-2.60 ± 1.02 deg (68% C.L.). This means that BICEP 
alone favors a non-zero Aa at about 2.4 a confidence 
level. Further, when WMAP7 and B03 data are added to 
the BICEP sample, the constraints get tightened again 
while the central value remains. In this case. 
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(2) 



FIG. 3: The one-dimensional posterior distributions of the 
rotation angle derived from the BICEP polarization data. 



This result gives a more than 3 a detection of a non- 
vanishing rotation angle. Note that we do not include the 
systematic errors of CMB measurements in our analysis, 
since it is not very clear how to combine those systematic 
errors together in a global analysis. 

Compared with WMAP7 and B03 data, BICEP data 
have smaller error bars, making it dominant in the joint 
analysis. Therefore, our result is largely due to the BI- 
CEP TB and EB polarization data. As an illustration, we 
plot a curve predicted by a Aa ~ —2.62 deg model with 
the data points in FiglT] Apparently, the bump struc- 
ture in BICEP TB data {£ < 400) is perfectly fitted by 
the curve, and another excellent fit can be found in the 
EB panel. In Figl3l we show the constraints on the ro- 
tation angle Aa from BICEP polarization data. We can 
see that the tight constraint on Aa mainly comes from 
the TB and EB power spectra of BICEP data, and TB 
and EB data give consistent limits on the rotation angle. 

Note, however, the sources of the CMB polarization, 
especially for the B-mode, are not unique. For example, 
the B-mode can be generated by the cosmological bire- 
fringence as mentioned above; it might be converted from 
E-mode by cosmic shear it could be the signature 
of the gravitational waves; or it can even be produced by 
the instrumental systematics . Therefore, one should 
bear in mind that the rotation angle might be degenerate 
with other cosmological parameters or nuisance parame- 
ters when fitted to the polarization data. As an example, 
we illustrate this degeneracy on the CMB BB power spec- 
tra in Fig|4l As shown, the three curves stand for three 
different mechanisms producing the BB power spectra. 
Interestingly, we find that the cosmological birefringence 
degenerate with tensor mode perturbations and with the 
gravitational lensing on large scales and on small scales, 
respectively. Therefore, in order to distinguish these ef- 
fects and obtain the clean information of the primordial 
tensor B-mode, the rotation angle has to be constrained. 
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FIG. 4: The theoretical predictions of the BB power spectra 
from three different sources: primordial tensor B-mode with 
r = 0.01 (black solid line); lensing-induced (red dashed line) 
and rotation-induced (blue dash-dot line). The cosmological 
parameters used here are Q,bh^ = 0.022, Qch^ = 0.12, r = 
0.084, Us ^1, As^ 2.3 X 10~^ and h = 0.70. 



and the measurements of TB and EB power spectra are 
really necessary. 

The QUaD data, shown in FiglU seem to have good 
quality, thus it is interesting to re-do the analysis on 
the rotation angle using this data set. We use the dif- 
ferent combinations of the two-point functions measured 
by QUaD and obtain the result shown in FiglSJ We start 
with the TE and EE data, and the constraint on Aa 
is found to be Aa = 0.01 ± 3.89 deg (68% C.L.). This 
constraint is very weak and the one-dimensional distri- 
bution is almost symmetric around Aa = 0, which is 
as expected - the rotated TE and EE spectra are re- 
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FIG. 5: The one-dimensional posterior distributions of the 
rotation angle derived from the QUaD polarization data. 



the time component of p^, will be much stronger. And in 
this paper we report a 3 ct detection of a non- vanishing 
rotation angle based on a joint analysis of the BICEP, 
B03 and WMAP7 data, when neglecting the systematic 
errors. 
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lated to the even functions of Aa as shown in Eq.(IT|), 
therefore the distribution must be invariant under a sign 
flip of Aa. Then we add BB data, the distribution be- 
comes bimodal, and peaks at |Aa| ^ 5.27 deg. The 
symmetry is also due to the foregoing arguments since 
BB still has no sign sensitivity on Aa. But the shift of 
the peaks is non-trivial. After a careful investigation, we 
find that the shift is due to the \ow-£ BB data, specially 
£ ^ 370, which might suffer from the unaccounted sys- 
tematic errors (23| . Then, we add all the two-point an- 
gular power spectra of the QUaD data together, and we 
find, Aa = 0.59±0.42 deg (68% C.L.), which is consistent 
with the published result in Ref. [l3l , but it seems in ten- 
sion with the results using other data samples, namely, 
QUaD is currently the only data sample favoring a pos- 
itive Aa at 68% confidence level. Finally, we combine 
QUaD with WMAP7, B03 and BICEP data, and the con- 
straint is further tightened to. 



Aa = -0.04 ± 0.35 deg (1 cr) 



(3) 



This means that the non-vanishing Aa preference dis- 
appears, which is because QUaD prefers a positive Aa, 
while other samples favor a negative one. 

Summary - In this work, we utilize the most recent ob- 
servational data of the CMB polarization to constrain the 
rotation angle Aa, an indicator of the cosmological CPT 
violation. Our results significantly improve the previous 
constraints on the rotation angle, and are already more 
stringent than those obtained from the polarization data 
of radio galaxies and quasars j24| . Furthermore we em- 
phasize that the radio galaxies and quasars only measure 
the rotation of the photon polarization from redshift up 
to z '-^ 2 till the present epoch, whereas for CMB polar- 
ization data the distance is much longer. Therefore, the 
constraint on the CPT violating parameter which is 
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